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ABSTRACT 

A number of factors that influence spectral position of the femtosecond pulse in a Kerr-lens modclocked CnLiSGaF 
laser have been identified: high-order dispersion, gain saturation, reabsorption from the ground state, and stimulated 
Raman scattering. Using the one-dimensional numerical model for the simulation of the laser cavity, the relative 
contributions of different factors have been compared. The Raman effect provides the largest self-frequency shift from 
the gain peak (up to 60 ran), followed by the gain saturation (~ 25 nm), while the high-order dispersion contribution 
is insignificant (~ 5 nm). Comparison with the experimental data confirm that the stimulated Raman scattering is 
a main cause of the ultrashort pulse self-frequency shift observed in Cr:LiSGaF and Cr:LiSAF lasers. 

Keywords: ultrashort pulses, Kerr-lens mode-locking, stimulated Raman scattering 

1. INTRODUCTION 

A rapid progress in Kerr-lens mode-locking technique allows to reach 14 and 12-fs pulse durations in Cr:LiSGaF and 
Cr:LiSAF lasers, respectively P^l These active media are attractive due to possibilities of sub-20 fs pulse generation 
directly from the diode-pumped lase r. At the same time, the lasers demonstrate a significant Stokes shift of the pulse 
spectrum at such short pulsewidths .EEE This shift decreases the accessible bandwidth due to the worse overlap of 
the gain and the pulse spectra, thus setting a limit to pulse duration .™ Therefore, investigation of the nature of the 
ultrashort pulse spectrum transformation has not only academic but also practical significance. 

A number of explanations for the ultrashort pulse spectrum shift in mode-locked lasers have been suggested. For 
example, it was supposed, that due to the high-order dispersions the spectral region of negative group-delay dispersion, 
which is optimal for pulse formation, may be displaced However in the framework of the perturbation theory 
the third-order dispersion does not contribute to the pulse carrier frequency, but results in spectrum asymmetry and 
fragmentation P Only in the case of the non-steady-state operation there is a possibility for the strong dynamic 
frequency shift P 

A more realistic explanation of the frequency shift takes into account frequency filtering due to reabsorption in 
the gain medium ? However, as it will be shown below, this explanation can not explain the dependence of the 
frequency shift on the pulse energy. Moreover, such dependence suggests that a nonlinear mechanism is involved in 
the frequency shift. 

As it was shown in Ref.P the stimulated Raman scattering in active medium is a very suitable nonlinear process, 
which can produce the experimentally observed Stokes frequency shift in sub-50 fs domain. The analysis in RefPwas 
based on the soliton perturbation theory that did not allow to take into account the effect of high-order dispersion 
and frequency-dependent dissipative lasing factors such as asymmetric gain band, spectral filtering on the output 
mirror, reabsorption in the gain medium and gain saturation. 

In this work we performed a numerical analysis of the spectral characteristics of ultrashort pulses on the basis of 
a relatively simple model of the Kerr-lens mode locking (KLM). The main advantages of our model is taking into 
consideration of the high-order dispersion, exact profiles of the loss and gain bands, frequency dependent reabsorption 
in the active medium, gain saturation and fast absorber action of the Kerr-lensing. The obtained results are in a 
good agreement with experimental data and allow to estimate the contribution of the different factors to spectral 
characteristics of ultrashort pulses. 
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Figure 1. General scheme of a KLM Cr:LiSGaF laser used in this paper. This scheme directly corresponds to the 
experiments in RefsPS HR, high reflector. CM, chirped mirror. OC, output coupler. 

The paper is organized as follows: first a summary of relevant measurements and experimental results is given, 
followed by the construction of the analytical and computational models. We then present the results of our calcu- 
lations and discuss the influence of each of the abovementioned factors separately. Finally, we present the results of 
simulation with all factors included, using distributed and discrete-element models. 

2. EXPERIMENTAL OBSERVATIONS 

Systematic Stokes shift of the u ltrashort modelocked pulse in Cr:LiSGaF- and Cr:LiSAF-based KLM oscillators has 
been first reported in 1997 GEEI where the stimulated Raman scattering has been suggested as a possible mechanism. 
It has been found, that i) the shift could be observed in oscillators with different dispersion characteristics, ii) the 
shift increases with the pulse energy, and iii) the shift increases with pulse shortening. At pulse durations below 20 
fs, the peak of the pulse spectrum may shift as far as 50-70 nm into the infrared as compared to the cw wavelength or 
modelocked spectrum at long pulse durations (^840 nm in Cr:LiSGaF) P Later, analogous shifts in Cr:LiSAF have 
been reported by Uemura and Torizuka^ and R. Gabel et al P^ All mentioned experiments used the common optical 
scheme, differing only in pump arrangements and dispersion compensation techniques. The schematic diagram of 
the laser oscillator is shown in Fig. ^ This is a representative scheme for any X-shaped KLM laser, because 
different types of dispersion compensation can always be represented by lumped dispersion of a chirped mirror. For 
modeling purposes we used the parameters of experiments, reported in RefsP^J the Brewster-angled 4 mm long 
LiSGaF crystal doped with 1.4% Cr, beam diameter in Cr:LiSGaF crystal 40x60 /.tm. The high reflectors (HR) had 
negligible dispersion. The dispersion of the chirped mirrors (CM) has been calculated from its original design, and 
additionally measured by white-light interferometry. The intracavity dispersion was calculated using the dispersion 
data of Cr:LiSGaF p measured dispersion curves of the chirped mirrors and calculated dispersion of the output 
coupler (Fig. [5J • We also used the experimental loss spectra (mirror transmission and ground-state absorption in 
CrrLiSGaF) as shown in Fig. 

The laser was pumped by 1.2-1.5 W from a Kr + -ion laser at 647 nm, with a TEMoo beam, of which 0.9-1.1 W 
have been absorbed in the active medium, generating 60-100 mW of average output power in the modelocked regime. 
The resonator was slightly asymmetric, with distances between the curved mirrors (radii of curvature 100 mm) 
and end mirrors being 88 and 109 cm, corresponding to 72 MHz pulse repetition rate. Modelocking has been 
achieved primarily by the hard aperture in form of an adjustable slit close to the output coupler (Fig. GJ. Fig. 01 
shows normalized output spectra at different pulse peak power, demonstrating the spectral shift. Similar spectral 
behaviour has been observed also in CrrLiSAF oscillators. However, for the sake of simplicity we provide experimental 
data and perform simulations on Cr:LiSGaF only. 

The Raman gain spectrum of undoped LiSGaF crystal has been measured according to the procedure described in 
Ref. P using the orientated crystalline quartz as a reference, and taking into account the thermal phonon population 
factor. The spectrum (see Fig. EJ) is obtained from spontaneous Raman scattering spectrum, recordered in confocal 
back-scattering geometry. Both incident and scattered light is polarized along the c axis, corresponding to the 
polarization of light in the laser. LiSGaF possesses spatial symmetry group (P31c) with 2 formula units in a unit 
cell, resulting in total of 32 optical phonon modes 3Ai g + 4:A 2g + kA\ u + 5A 2u + 8E g + 8E U , of which 3Ai g + 8E g are 
Raman-active. In the scattering geometry as described above, only 3 full-symmetric Ai g modes are visible (Fig. EJ, 
with relevant parameters given in Tabled 
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Figure 2. Measured group delay dispersion of the active media (8 mm in double pass), output coupler and chirped 
mirrors in dependence on the wavelength A. 

Table 1. Raman gain of undoped LiSGaF. 



Frequency f2j/2-7r 
cm -1 


Raman gain 
cm/GW 


Width (FWHM) 
cm -1 


T 3 
ps 


230 


0.014 ±0.005 


9±3 


1.2 ±0.4 


349 


0.021 ±0.006 


14 ±2 


0.7 ±0.1 


551 


0.32 ±0.05 


12.5 ±0.6 


0.86 ±0.05 



3. MODEL 

There exist different approaches to modeling of ultrashort pulse generation in solid-state laser, which are based 
on soliton or fluctuation models. The soliton approach can be applied only the distributed laser model (where 
the dispersion and nonlinearity are implied to be evenly distributed over the round-trip and act simultaneously) 
but allows to build comparatively simple analytical description thus promoting easy interpretation of results. We 
also based our calculations on the distributed laser model but the results were tested by simulations on the basis 
of discrete-element scheme corresponding to Fig. Q] To overcome the limitations of the soliton approach we used 
numerical simulations allowing to account for the high-order dispersion, the laser field reabsorption, the complicated 
spectral profiles of the gain and output coupler transmission, and the Raman scattering within the active medium. 

The modelocking is described by a fast absorber-like action of Kerr-lensing in the active medium in the form of 
a nonlinear transmission-operator exp — 1+tJ ^( z t)\ ,j! ' where 7 is the modulation depth (KLM loss), which is set 
by the cavity arrangement, a is the inverse loss saturation intensity, a is the field depending on local time t and 
longitudinal coordinate z (\a\ 2 has the meaning of the field intensity). Parameters 7 and a are controlled by changing 
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Figure 3. Ground-state absorption of Cr:LiSGaF (a), round-trip resonator losses due to the output coupler (b) 
and gain cross-section of Cr:LiSGaF (c). 

the cavity configuration, which is a common procedure for Kerr-lens mode-locking experiment. 

Another fundamental factor in our model is the presence of high-order dispersions due to the active medium, the 
dispersion compensator, the output coupler, and high-reflective mirrors. The corresponding experimental character- 
istics are shown in Fig. [5J For the numerical calculations, the data were represented by the eighth-order polynomial 
approximation. The action of dispersion can thus be presented in the following form: 



a(z,t) = / a(z,t )G(t - t')dt ', 



(1) 



G(t-t') = ^- J exp l-iJ2^ D j(v-uoy - (t-t')ujjdoj, 



(2) 



where ui is the frequency, G(t — t') is the Green's function depending on the dispersion coefficients Dj up to eighth 
order of j. 

The next important factor in our model is the gain saturation that was described on the basis of quasi-two level 
scheme of the active medium operational 




Figure 4. Experimental demonstration of the pulse frequency shift. The output spectra of modelockcd Cr:LiSGaF 
laser, with changing intracavity pulse peak power (bottom spectrum corresponds to the lowest peak power, top 
spectrum - to the highest peak power). 
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Figure 5. Raman gain of undoped LiSGaF. Exciting laser line at 514.5 nm and scattered light are both polarized 
along the crystallographic z axis, corresponding to the polarization of the laser radiation in Cr:LiSGaF laser. 



da IpcriA , . \a\ 2 cr 32 a 

where I p is the pump intensity, v is the pump frequency, CT14 and (T32 are the loss and the gain cross-sections, 
respectively, T r — 85 /is is the gain relaxation time. If the pulse duration is much less then cavity round-trip time 
T cav = 14 ns, then this equation can be replaced by the following one: 

- = P {amax - a) -- a -— a (4) 

where P — {I p <Ju/hv)T cav is the dimcnsionless pump intensity, z is the dimensionless longitudinal coordinate, i. e. 
the number of the cavity round-trips. E s — hv/a^2 is the gain saturation energy flux, E is the full pulse energy flux. 

Self-phase modulation in the active medium was represented by nonlinear "transmission" operator exp —i(3 \a(z, t)\ 2 

where (3 = 2irri2x/\n — 3.4 cm 2 /TW is the self-phase modulation coefficient. Here n and 77,2 are the linear and non- 
linear coefficients of refraction, respectively, A is the central wavelength corresponding in our case to the gain band 
maximum, x = 8 mm is twice the length of the active crystal. 

Finally, we consider the stimulated Raman scattering within the active medium. Following Ref. P where Raman 
scattering contribution was calculated analytically on the basis of the soliton model, we supplement the model with 
the following equations: 



da s 
da p 



3 



= i ^2Q* 1 a p> ( 5 ) 

i=i 

3 
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dt 2 Ti dt 



fljQj = /JjOpO*, (7) 



where <; is the longitudinal coordinate (pulse propagation axis) inside the active medium, a PyS are the amplitudes 
of the "pump" and the "Stokes" components within generation spectrum, fij are the phonon resonance frequencies 
(j = 1,2,3 corresponding to the three Raman-active phonon resonances in LiSGaF, see Fig.EJ). Tj are the inverse 
bandwidths of Raman lines, fij—g^flj /Tj are the coupling parameters for Raman gain coefficients g?. 

Solving third equation of the system, we obtain the steady-state phonon amplitude for the fixed pump and the 
Stokes components with corresponding frequencies lu p and uj s : 



p. [ija p a s ^ (ija p a s . . 

ttj T] ~ [Up - u s ) 2S2j (Uj - (lo p - Lu s )) - - T ^- 

The validity of the last approximate expression follows from the fact that the Raman lines are narrow in compar- 
ison with the pulse spectrum. With regard to the contribution of all spectral components of the pulse to the phonon 
amplitude, the equation for the Stokes and the pump fields in the frequency domain can be written as: 

3 2 
oa s _ \ ■* \ ■* \ a p,k\ ,„s 

° q 3=1 k 2 "j(H/ _ ( W M -U s ))+ -jT 1 - 

3 2 

da P _ \ -» \ -> \a s ,k\ / ln s 

j=1 k 2n J (n j -(ujp~uj s , k ))-^ F f- 



where k is the index of the field's spectral (i.e. Fourier) component (in the simulations we considered 2 13 components). 
Since the Raman lines are narrow, the field variation within these lines is negligible, i.e. a SjP are constant and can 
be taken out of the second summation. Then the summation can be executed explicitly by transition to the integral, 
resulting in 



da s 

!k 4 
da p 

Ik = ~4 



J]T> ffj s |a p | 2 , (11) 



3 



= -iI>^H 2 . (12) 

3=1 



It should be noted, that in these equations u> p — lo s = Qj and there is the pair-wise interaction of the spectral 
components within the wide enough generation spectrum. 

There are two main mechanisms of the generation of the initial seed at the Stokes frequency for Eqs. The 
first one is the spontaneous Raman scattering with the increments of the scattered intensity growth 



Xj 



exp 



/ Mij y 

\~kFTj 



(13) 



where n s is the index of refraction at Stokes frequency, T is the temperature, and fcs is the Boltzmann's constant 
(HI With this seed signal, the stimulated Raman scattering results in appearance and growth of spectral replicas of 
the main oscillation pulse, shifted to the lower frequencies by the Raman frequencies Qj. 

More significant source for Stokes component's amplification, however, is the broad-band pulse field itself. When 
the pulse pulse spectrum is wide enough to become comparable with the Raman frequency shift, the lower-frequency 
part of the spectrum can play a role of the Stokes component seed with respect to the higher-frequency part of the 
spectrum. The stimulated Raman scattering then transfers the energy form the higher-frequency components to the 
lower- frequency ones, resulting in the continuous red-shift of the pulse spectrum as a whole. As the field amplitude 
of the laser pulse significantly exceeds the spontaneous seed, the second mechanism strongly dominates over the first 
one. However, we included both mechanisms in our model, because their influence on the pulse spectrum is quite 
different. 

Later on it is convenient to normalize the time to the inverse gain bandwidth t g = 2.25 fs and the intensity to 
resulting in the normalization of the field energy to (/3i s ) _1 . As already pointed out, we analyzed the described 
above model in two ways: on the basis of distributed and discrete-element approaches. In the framework of the 
distributed model, we didn't consider the propagation through the individual laser element and supposed that the 
pulse envelope is formed by the overall net-dispersion in the cavity. As result we have a split-step scheme describing 
ultrashort pulse propagation from z to z + 1 transits: 



a(z' 1 1) = 
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(14) 
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n(.:-r)=n(:l,.x 1 ,|-r / „l .:. t')f dt' - T cav /T r - P ) + ^ WT r P» (lg) 

* ~T J- cav / J-r 



where r = t g /(E s /3) = 0.00079 is the normalized gain saturation parameter. The Green's functions A, L, C describe 
spectral bands of gain, reabsorption and output loss, respectively (note, that the dispersion is already is included in 
G): 



A{t-t') = -!- / (l + a(z))^ a (u))exp[ioj(t-t')]du, (17) 
L(t-t') = / $ r (u) exp [ioj (t - 1')] du, (18) 

C(t - 1') = — j <I>,„„ U'icxp ['-'(' ■ I'f <U. CUD 

Here, r _ ou t are the "form-factors" describing spectral profiles of the gain, reabsorption and the output coupler 
bands, which resulted from the functional approximation of the experimental data (Fig. [2J ■ 

The system (0 has to be completed by the system (|2J and the result can be solved on the basis of numerical 
simulation in Fourier domain and split-step method for nonlinear propagation. 

The discrete-element approach is based on the element-to-element simulation of the pulse propagation through 
the cavity on every round-trip, following the laser scheme in Fig. ^ The nonlinear and spectral characteristics of 
each laser element are considered separately. Further refinement of the simulation is achieved by considering the 
pulse propagation through active medium by splitting it into five slices and applying the split-step procedure to each 
slice consecutively. 



4. DISCUSSION 

Our simulations are aimed at investigation of the influence of the different factors on spectral characteristics of the 
ultrashort pulses. Therefore, to simplify the interpretation, we will first consider the high-order dispersion action 
without Raman scattering and without reabsorption; then reabsorption and Raman scattering will be taken into 
consideration without high-order dispersion action; and finally, the join action of all factors will be analyzed. To 
conclude, we will also compare the obtained results with the calculation based on the discreet-element model. 

4.1. High-order dispersion action 

As it was discussed in p small contribution of third-order dispersion to characteristics of Schrodinger soliton does not 
cause the frequency shift, but introduces additional field time delay on the cavity round-trip. However, the influence 
of the linear and nonlinear dissipative terms in Eqs. J3J can destroy the soliton character of the pulse, in particular, 
it can add the chirp. The latter, as it was shown in P results in the frequency shift of the pulse spectrum in the 
non-steady-state regime. Additionally, in the real-world laser systems the contribution of the high-order dispersions, 
as a rule, lies beyond the bounds of the perturbation theory validity. 

The typical net-dispersion curves corresponding to the experiment with the chirped mirrors are presented in Fig. El 
a. Over the full spectral region of generation the pulse undergoes the non-negligible influence of the dispersion up 
to the eighth-order (see Eq. (J2J). As a rule, there is the closed spectral window of dispersion, which is "optimal" for 
steady-state pulse generation. The contribution of the high-order dispersion terms can shift this window. In Fig. [|J1 
a, this shift corresponds to the red-shift of the positive net-dispersion branches (transition from solid to dash and dot 
curves). The resulting output spectra obtained from the distributed model are shown in Fig. El b. The net-dispersion 
shift is accompanied by the red-shift of the pulse spectrum (transition from solid to dash and dot curves). Besides 
this effect there is the possibility of the essential spectral profile distortion (solid curve) and side-band generation 
(side-bands lie outside the of shown region) El 

However, we cannot consider this shift as the cause of the experimentally observed effect because the dispersion 
shift has the linear nature, i. e. there is no obvious dependence of this shift on the field energy. This is demonstrated 
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Figure 6. The dependence of dispersion (a) and generation spectra (6) on wavelength. P = 3.2x10 4 , a = 1, 7 
= 0.05, pulse energy E is 20 nJ. Pulse durations t p : 27 (solid curve), 38 (dash), 36 fs (dot). 

by Fig. |3 where the pulse energy variation due to the pump variation changes the spectral profile, but does not 
cause any noticeable frequency shift (compare solid and dash curves in this figure). 

The obtained results demonstrate that the self- frequency shift cannot be caused by the non-dissipative factors. 
As the pulse duration is too large for the nonlinear dispersion to play any significant role, we will concentrate on 
spectrally-dependent losses and Raman effect. 

4.2. Output loss and reabsorption in gain medium 

As mentioned above, dependence of the frequency shift on pulse energy implies the involvement of some nonlinear 
mechanism. Since in the real-world systems the gain band does not coincide with filtering band (output mirror in 
our case) and the reabsorption band, the spectral position of the net-gain maximum changes with the gain value. 
The latter is defined by the pump and by the pulse energy (see Eqs. iJSJ): pulses with higher energy experience 
lower gain due to the multi-pass saturation. The dependence of the net-gain maximum on the saturated steady-state 
gain coefficient a is shown in Fig. |H1 by the solid curve. This curve was obtained from the numerical analysis of 
the measured spectrum of the output coupler, intracavity loss and gain profiles. As we can see from this figure, the 
behavior of the net-gain maximum corresponds to the Stokes frequency shift with the pulse energy growth, due to the 
gain coefficient decrease. However, the magnitude of this shift (~ 10 nm) is not sufficient to explain the experimental 
values (up to 50 nm). Besides that, the location of the pulse spectrum in general does not coincide with the net-gain 
maximum. 
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Figure 7. The dependence of generation spectra on wavelength in the presence of high-order dispersion in the 
distributed (solid and dash curves) and the discrete-element (dot) models, a = 1, 7 = 0.05. For the corresponding 
pulse parameters see Table |21 




Figure 8. Pulse central wavelength as a function of the saturated gain coefficient. Solid line: Net-gain maximum; 
ABCD spectrum peak in the case of gain saturation without reabsorption and Raman scattering in the active 
medium; EFGH contribution of the reabsorption; UK contribution of Raman scattering. Points correspond to 
Table El 



Table 2. Normalized pump power, pulse duration and energy for Fig.[7| 
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dot 
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Table 3. Ultrashort pulse energies and durations for Fig. (HI 
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Energy, nJ 


16 


27 


47 


54 


20 


38 
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33 


78 


103 


Pulse width, fs 


45 


28 


20 


23 


28 


20 


18 


37 


25 


14 
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To check the last thesis we performed numerical simulations on the basis of the described above model in the 
presence of dissipative factors only. The results are presented in Fig.[S]by the dashed lines ABCD and EFGH. The 
letters alongside of squares, circles and triangles denote the intracavity pulse energy obtained from the numerical 
simulation (see Table 0J. The pulse energy variation results from the change of modulation parameter a, i.e. due to 
the change of the cavity configuration. 

The curve ABCD illustrates the case when there is no reabsorption and Raman scattering in the active medium. 
The dependence of frequency shift on pulse energy is not monotonous: at the minimum pulse width (15 fs in our 
case) there is no Stokes frequency shift, contradicting the experimental data. Moreover, the maximal shift (of about 
6 nm) is by an order of magnitude less than that observed in the experiment. It is clear, that gain band asymmetry 
and output loss alone cannot produce large self-frequency shift in Cr:LiSGaF laser. 

Stronger contribution is provided by the reabsorption in the active medium (curve EFGH on Fig.[5J). For small 
energies the simulation result (point E on Fig. |SJ) coincides with the prediction made on the basis of elementary 
model of net-gain maximum shift (solid curve) . The rise of the pulse energy increases the shift (up to 25 nm in our 
case), much more than in the previous case, although still by a factor of two less than in the experiment. 

Note, that after certain pulse energy the spectral shift is again decreasing. This is caused by the increase of the 
pulse duration for the large energies (point H on Fig. [SJ) , due to the nonlinear loss saturation, resulting in reduced 
spectrum width. The laser approaches the condition of the cw operation, described by solid curve (although the 
conditions of pulsed gain saturation strongly differ from those in cw-regime). 

4.3. Raman scattering 

As pointed out in the previous section, the net-gain shift model fails to provide correct description of the experimental 
data by a factor of two. However, taking the stimulated Raman scattering into account allows to obtain large 
frequency shifts, increasing with pulse energy growth. In Fig. [SI the curve UK demonstrates the Raman scattering 
action in the absence of reabsorption effect. As the simulation demonstrates, the red components originate from the 
amplified Raman signal, which pulls the whole spectrum over the long wavelength limit at the given pump power, 
defined by the spectral filtering. 

Note the pronounced threshold- like character of the effect. For the small energies (point /) the spectral shift is 
negligible, but the energy growth causes very strong shift (60 nm in our case) in good agreement with experimental 
results EEE1 anc i w ith analytical prediction W Since the gain saturation does not play important role in this case, 
the dependence of frequency shift on gain coefficient is insignificant. 

As already mentioned above, the main contribution to the stimulated Raman process comes from the energy 
transfer from the blue part of the pulse spectrum (pump) to the red one (Stokes) . The efficiency of the stimulated 
Raman scattering is therefore defined by the product of intensities at pump and Stokes frequencies. The separation 
between Stokes and pump components is fixed, it is equal to the Raman line frequency ttj. Therefore, decreasing 
the pulse spectrum width strongly suppresses the effect and reduces self- frequency shift. Assuming that the pulse 
spectrum has exponential fall-off to the blue and red sides, we see that the dependence of the Raman shift on the 



pulse spectrum width should be asymptotically exponential at long pulse durations. This is also the result of the 
analytical theory in RefP where sech 2 pulse shape has been assumed. 

As the Raman effect strongly depends on the pulse intensity and Raman gain, it should be especially pro- 
nounced in low-gain lasers working with low output coupling and high intracavity pulse energy, such as Cr:LiSGaF, 
CnLiSAF, Cr:YAG. These materials also possess stron g and broad Raman lines P Large power-dependent red-shift 
in Cr:LiSGaF and Cr:LiSAF is well documented EE0 l n femtosecond Cr:YAG lasers, femtosecond pulse spectrum 
is also always red-shifted with respect to the cw wavelength in the same resonator 

4.4. Discrete-element model 

Finally, we can compare the simulation results in case of distributed and discrete-element models. In the case of 
the high-order dispersion action the transition to the distributed-element model does not significantly change the 
spectral characteristics of the pulse (the dot curve in Fig. El b). The long- wavelength "shoulder" of the spectrum in 
the case of the net dispersion corresponding to solid curve in Fig. [fj] a is somewhat stronger than in the distributed 
model. This is because this "shoulder" results from local dispersion maximum due to the chirped mirror (CMi on 
Fig. 0, which is the closest to the output mirror. Additionally, the self-phase modulation contribution is found to 
be slightly higher in comparison to the distributed model. However, all these changes are not qualitative. We found 
also that the transformation of the spectrum on each laser elements is small ensuring the validity of the distributed 
model. It should also be noted, that the transition from the distributed to the discrete model slightly increases the 
contribution of the Raman scattering, seen by the lower threshold energy. 

5. CONCLUSION 

Using the numerical simulations performed in the framework of one-dimensional distributed and discrete-element 
models, we analysed the spectral characteristics of a cw Kerr-lens mode-locked Cr: LiSGaF-laser. The two main 
factors causing the ultrashort pulse self-frequency shift have been established: the nonlinear shift of the net-gain 
band due to the gain saturation in the presence of reabsorption in the active medium and the stimulated Raman 
scattering. The first effect is essential for comparatively small pulse energies and produces wavelength shifts up to 30 
nm. The Raman scattering occurs as a result of pulse energy growth and causes the large (over 50 nm) red shift . The 
contribution of the high-order dispersion (up to the eighth order), gain-band asymmetry and spectral characteristics 
of output coupler were estimated as well. However, their effect on the pulse central frequency is much smaller than 
that of the stimulated Raman scattering, which is the main cause of spectral red-shift in Kerr-lens modelocked laser. 
The shift values obtained from the numerical simulations are in good agreement with experimental data. 

Analytical part of the calculations is presented on http://www.geocities.com/optomaplev 
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